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Abstract 



The renormalisation group running of the quark mass is determined non-perturbatively for 
a large range of scales, by computing the step scaling function in the Schrodinger Func- 
, tional formalism of quenched lattice QCD both with and without 0(a) improvement. A 

one-loop perturbative calculation of the discretisation effects has been carried out for both 
the Wilson and the Clover- improved actions and for a large number of lattice resolutions. 
The non-perturbative computation yields continuum results which are regularisation inde- 
ed I pendent, thus providing convincing evidence for the uniqueness of the continuum limit. As 
a byproduct, the ratio of the renormalisation group invariant quark mass to the quark mass, 
renormalised at a hadronic scale, is obtained with very high accuracy. 



1 Introduction 



The renormalisation group running of the QCD fundamental parameters, namely 
the renormalised gauge coupling and quark masses, has now been computed non- 
perturbatively for a large range of scales, albeit in the limit of infinitely heavy sea 
quarks; see refs. These results have been obtained using lattice regularised 

quenched QCD with Wilson fermions, prior to extrapolating to the continuum limit. 
In the case of the quark mass running, Symanzik improvement was an important 
element in reducing the extrapolation uncertainties, as it implies that the dominant 
systematic effects due to the finiteness of the UV cutoff are O(a^). The case of 
QCD with two dynamical flavours has been investigated in '3 . 

The continuum limit of lattice QCD is known to exist to all orders of pertur- 
bation theory (PT) 4 . Beyond PT this issue has been addressed by numerical 
simulation. The strategy consists in implementing different regularisations which 
formally correspond to the same Field Theory (QCD in our case) and in establish- 
ing the universality of the continuum limit of given renormalised physical quantities, 
computed with different regulators. This sometimes turns out to be less straight- 
forward than expected; e.g. see the discussion on the universality of the continuum 
limit of spin and sigma models in refs. In pure SU(2) gauge theory, universal- 
ity has been tested by computing two different non-perturbatively defined running 
couplings over a large range of energies |n|. In pure SU(3) gauge theory, evidence 
of universality has been recently found in a study of the scaling properties of the 
deconfining temperature with different gauge actions 0. 

In the present work we extend these ideas to the step scaling function (SSF) 
of the quark mass in quenched QCD. This quantity has been calculated in ref. [2] 
from an improved action; here it is also evaluated from an unimproved action. The 
continuum SSF, computed for a large range of renormalisation scales, is found to 
be independent of these regularisation details, providing evidence for a universal 
continuum limit. In this respect our study parallels closely the one of ref. [S], which 
was dealing with the SSF of the operator corresponding to the average momentum 
of non-singlet parton densities. In comparison, our quantity is particularly simple, 
as we essentially compute ratios of two-point functions. This allows us to have an 
excellent control of both statistical and systematic errors. As a byproduct we recal- 
culate the flavour independent ratio of the renormalisation group invariant (RGI) 
quark mass to the renormalised one (at a given hadronic scale). We obtain a result 
compatible to the original one of ref. but fairly more accurate. 

A study of discretisation effects in the SSF has also been performed, for both 
Wilson and Clover actions, in one-loop PT. This calculation has been carried out for 
a large number of lattice resolutions. We find that lowest order perturbation theory 
greatly underestimates the discretisation effects of the SSF. 
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2 The Schrodinger Functional and 0(a) improvement 



In this section we gather the most relevant definitions and outhne the properties 
of the quantities we are interested in. Most details are omitted, as they have been 
presented in previous works, which we will frequently refer to. 

We adopt the lattice Schrodinger functional (SF) formalism |9| lin | ITT ] : more 
specifically we regularise QCD on a lattice of extension x T (here T = L al- 
ways) with periodic boundary conditions in the space directions (up to a phase 9 
for the fermion fields) and Dirichlet boundary conditions in the Euclidean time di- 
rection jlUl lllj . Otherwise the lattice gauge and fermionic field actions are of the 
standard Wilson type; their 0(a) improved version is discussed below. The bare 
gauge coupling and quark mass are denoted by go (with (3 = Q/qq) and mo (with 
2k = [arriQ -\- 4]~^), respectively. As we will be working in the quenched approxi- 
mation, the bare gauge coupling go and chiral point Kc are functions of the lattice 
spacing a alone. The chiral point is the value of the hopping parameter k for which 
the "current" quark mass, defined below, vanishes. The bare subtracted quark mass 
is defined as aniq = [1/k — 1/kc]/2, whereas an unrenormalised "current" quark 
mass is given by 

"^^^"^ = 2/p(xo) ' ^'-'^ 
with /x (X = A, P) the correlation functions of local bilinear operators 

t 
2 



/x(xo) = 5^(0(y)75Ci(z) V^i(x)7xV,(x)) . (2.2) 



The field indices i,j label two distinct flavours; the "boundary fields" ( are defined in 
ref. ^Ij. For X = A we have 7x = 7o75 and for X = P we have 7x = 75- The forward 
and backward lattice time derivatives are denoted by do and 9g respectively^. We 
also define the correlation function of boundary fields 



/i = -^ E (c;(y')75Ci(z')0(y)750(z)) . (2.3) 

y.z.y'.z' 

Unprimed quantities are defined on the the xq = boundary, primed ones on the 
xq = T one. 

The 0(a) Symanzik improvement of the above construction has been worked out 
in refs. |in|ll2j. For the pure gauge action, it amounts to modifying it by introducing 
time-boundary counterterms proportional to [ct(go) ~ Foi^ the fermionic action we 
must introduce the well-known clover counterterm in the lattice bulk, proportional 
to Csw(fl'o)) time-boundary counterterms proportional to [ct(g'o) — 1]. Correlation 

^We follow closely the notation of ^1 112|. whither we refer for details. 
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functions of composite operators such as eq. (|2.2|) may then also be 0(a) improved 
by including in the lattice definition of these operators the appropriate counterterms. 
In the chiral limit there are no such counterterms for the pseudoscalar density P{x), 
while the axial current Ao{x) requires the addition of doP with a coefficient Ca(5o)- 
The axial current is used in the computation of the bare quark mass, but, being 
scale independent, it is clearly not needed in the computation of its renormalisation 
group running. Thus Ca will play no role in the present work. 

All these improvement coefficients may in principle be computed non-perturba- 
tively for a range of values of the bare coupling go; for Cgw we rely on the calculation 
of ref . ^31 • It has also been calculated in perturbation theory to one loop jl4| II 5j . 
The coefficients Ct and Ct are known only in perturbation theory, to NLO [16j and 
LO ^H] respectively: 



A more detailed discussion of perturbative 0(a) improvement will be presented 
in Sectional Here we outline the main expectations related to cutoff effects, in the 
spirit of the Symanzik improvement programme ^\ : in the absence of improvement 
counterterms (i.e. Cgw = and ct = Ct = 1), correlation functions (such as fp and 
/a), computed at fixed UV cutoff and renormalised non-perturbatively, should 
exhibit 0(a) deviations from their continuum limit. If all improvement coefficients 
were known non-perturbatively, the discretisation errors would be 0(a2). With the 
improvement coefficients set to their tree- level values (i.e. Csw = Ct = Ct = 1), the 
dominant discretisation effects are expected to be 0(^0^) and O(a^); with one-loop 
coefficients we have 0{gQa) and O(a^) errors etc. These statements refer to the 
chiral limit (away from which, we must also take into consideration counterterms 
proportional to the quark mass). Since we are in the framework of mass independent 
renormalisation, working in the chiral limit is adequate for our purposes. We have 
performed numerical simulations in two regimes: 

(i) What we call "unimproved action results" (or "unimproved case" for short) 
consists in setting Csw = . Moreover, we set ct = 1, while the one-loop value^ 
(eq. (|2.4|) truncated to 0(go)) is used for ct- Since the action in the lattice 
bulk is unimproved, the dominant discretisation effects ought to be 0(a). 

(ii) What we call "improved action results" (or "improved case" for short) consists 
in using the Clover action with a non-perturbative Cgw The one-loop value 

^This is a choice of convenience: it is important to know for renormalisation purposes (see 
eq. below) the dependence of the Schrodinger functional renormalised coupling g{l/L) on 

the bare coupling go- This dependence is known non-perturbatively for the pure Yang-Mills 

action with this Ct value. In any case, the choice for Ct has no bearing on the order of leading lattice 
artifacts. 
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from eqs. (|2.4|) and (|2.5|) is used for ct and ct respectively. Thus the domi- 
nant discretisation errors should be 0{gQa) and O(a^). Since the former only 
arise from perturbatively improved boundary counterterms, while everything 
in the lattice bulk is fully improved, it is reasonable to expect that correla- 
tion functions are mostly affected by O(a^) errors. Numerical support for this 
expectation has been presented in ref. 



3 The step scaling function 



SF renormalisation schemes are mass independent; i.e. simulations can be performed 
in the chiral limit. The renormalisation scale is set at the lattice IR cutoff (i.e. 
fi = l/L); the renormalised coupling g{\/L) and quark mass m{\/L) are then only 
functions of L. The SF renormalised coupling has been defined in ref. pp. The 
renormalised quark mass is [2] 



m{l/L) = lini Zj^{gQ) Zp^{go,L/a) m{go) 



(3.1) 



where m{gQ) is defined in eq. (|2.1() and the renormalisation condition for the pseu- 
doscalar operator is 



Zp{go,L/a) 



ML/2) 



c{9,a/L) , 



(3.2) 



m=0 



with c{9, a/L) such that at tree level Zp(0, L/a) = 1. We will always impose eq. (|3.2|) 
at 6 = 0.5 [21 and hence eliminate any explicit reference to 9 from now on. 
The axial current normalisation Z^{gQ), being scale independent, has no effect on 
the renormalisation group running of the quark mass; thus it is of no immediate 
consequence to the present work. 

Here we are interested in the step scaling function of the quark mass, which is 
defined in the chiral limit rra(go) = 0, for a lattice of a given resolution L/a and at 
fixed renormalised coupling g'^{l/L) = u, by 



Sp(u, a/L) 



Zp{go,2L/a) 



Zp{go,L/a) 

This quantity is finite in the continuum limit 



(3.3) 



m=0, g^{l/L)=u 



/ ^ 1- V / /T\ m{l/L) 
dp u = liniSp(u, a/L = _ 

a^o m[l/{lL)) 



(3.4) 
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The physical meaning of dp follows from the RG equation obeyed by the renormalised 
quark mass 



t^^-^ = r{m) (3.5) 

(recall that /i = 1/-L). Upon integration of this equation between scales and 
(2L)^^ we obtain 

ap(u) = exp <^ - / — - dg ) , (3.6) 
[ Jg{l/L) P[9) J 

with P{g) the Callan-Symanzik function. Thus, ap is closely related to the quark 
mass anomalous dimension. 

The lattice SSF Sp is not unique: it depends on the details of the lattice 
regularisation (e.g. the type of lattice action chosen, the level of 0(a) improvement 
etc.). Its continuum limit, however, should be unique (i.e. universality should hold), 
unless lattice QCD, or at least the specific regularisation implemented here, exhibits 
some unexpected pathology. This is what the present paper has set out to explore, 
in the spirit of refs. [SlEl- 



4 Discretisation effects in perturbation theory 

The expansion of the SSF in renormalised perturbation theory reads 

oo 

ap(7x) = l + ^4")^x" , (4.1) 



n=l 



with the LO universal RG coefficient = — 81n(2)/(47r)^. In perturbation theory 
the cutoff dependence of the SSF can be studied by expanding 

M^ = l + f;fc„(a/L)^V (4.2) 

Note that due to the choice of renormalisation condition ()3.2() . discretisation errors 
are absent at tree level (i.e. kQ^a/L) = 0). Moreover, kn{0) = by construction. 
The quantity ki{a/L), which contains the cutoff effects at one loop, is known for the 
improved case from the work of Sint and Weisz |18j for L/a = 4, 6, . . . , 16. In their 
notation, it is given by 



ki{a/L) = k{^)6k{L/a) = -^^^6k{L/a) , (4.3) 
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L/a ki (Unimproved) ki (Improved) 
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Table 1 : Results for the cutoff dependence of the step scaling function of the pseu- 
doscalar density at one loop in perturbation theory. 

with 6k{L/a) tabulated in Table 2 of ref. |1^ (the case of interest to us is ^ = 
0.5, p = T/L = 1). We have repeated these calculations for the unimproved case 
and extended both cases to L/a = 18,20,... ,26. The results are summarised in 
Table Hand Fig.IH 

The present perturbative analysis has been motivated by the wish to explore 
in detail the cutoff dependence of our non-perturbative estimates of Ep(?/, a/L), 
obtained in our simulations at four lattice resolutions L/a = 6,8,12,16 (see next 
Section). Strictly speaking we therefore need to know ki{a/L) only at these four 
values of {a/L), for both the unimproved and improved cases. However it is clear 
that one- loop cutoff effects can have a rather non-trivial overall behaviour as (a/L) 
is reduced. Had we limited our perturbative calculation to the range of interest (i.e. 
L/a = 6,8,12,16), in the unimproved case we would have only observed that ki 
crosses over monotonically the abscissa axis at about L/a = 8, without signalling 
that it indeed converges towards its limiting value A;i(0) = (the improved case 
already "bends over" towards the origin within this range). As a partial safeguard 
against the eventuality of some uncontrolled error afflicting our perturbative results 
(e.g. rounding in the numerical integrations), we have extended the calculations all 
the way to L/a = 26. Tableland Fig. ^ demonstrate that indeed ki reaches a local 
extremum and subsequently points towards the origin of the it should. 

In order to gain some further insight into this behaviour, we recall that the 
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coefficients fc„ of tlie perturbative series can be expanded as 



K{a/L) = {jY T.'^Pi^ IH'^/L)]' , (4.4) 

p=l 1=0 

with the leading 0(n) (one-loop) discretisation effects having the form 



hia/L) 



L 



a" 



+ cS? ln(a/L) 
+ ln(a/L) 



(4.5) 



Tree-level improvement implies in general that c^^'' = 0; thus the 0(a/L) one-loop 
perturbative contribution is as above, but without the logarithm (i.e. c^^ = 0). 



One-loop improvement implies that c^"'' = c 



l,n-l 



0; thus the 0{a/L) one-loop 
0). The leading 



(1) 

10 



"11 



perturbative contributions of eq. (|4.5|) all vanish (c 
discretisation effects of Sp are then ©(ua^/L^). Such "dominant" discretisation 
effects may in practice compete with the next order 0(u^a/L) errors, arising from 
k2{a/L). In all cases the functional form of ki{a/L) is clearly complicated; the 
observed behaviour can be explained by the strengths (and relative signs) of the 



various coefficients c, 



(n) 
pi 



5 Non-perturbative computation of the step scaling function 

In this Section we study the extrapolation of Sp to the continuum limit. We also 
obtain a very accurate estimate of the ratio of the RGI quark mass to its renormalised 
counterpart at a hadronic scale. The method of computation is identical to that of 
ref. [2]. 

5.1 Continuum limit of the step scaling function 

For both the unimproved and improved cases the lattice SSF Sp has been evaluated 
at 14 values of the renormalised coupling g{l/L), each for four lattice resolutions 
L/a = 6,8,12 and 16. Note that in ref. |2] only the improved case has been studied. 
A full collection of our raw data is presented in Tables |21 and IHI The "tuning" of 
P at the four L/a values, corresponding to an (almost) fixed renormalised coupling 
5^(1/L) = u has been taken over from ref. 0. This corresponds to the first three 
columns of Tables [^J and |21 The same is true of Hc for the improved case (fourth 
column of Table IS}. All other results are new^. In the strong coupling regime new 

^The computation of Kc for the unimproved case was first performed in ref. [H| for 9 of the 14 
couphngs used here; following this is done at 9 = 0. 
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Figure 1: Cutoff dependence of the step scaling function of the pseudoscalar density 
at one loop in perturbation theory. 

results have been obtained at g'^{l/L) = 3.111. The statistical accuracy of our 
improved and unimproved Sp(u,a/L) results is comparable (save for a few cases 
where the improved data has somewhat smaller errors).^ 

A comparison of our data in the improved case with those of ref. reveals 
fairly compatible results: for Zp{gQ, L/a) and Zp{gQ,2L/a) we mostly agree within 
errors, save for a few cases in which agreement is within 2a; the same is true for 
Sp(u,a/L), with some data being compatible only within 1.5 o". 

If quenched lattice QCD has a universal continuum limit, then both sets of 
Sp results (improved and unimproved action) ought to extrapolate to the same 
continuum value dp at fixed coupling u. What we have set out to investigate is 
the power dependence (linear and/or quadratic) of the results on {a/L). From the 
analysis of ref. we expect the dominant discretisation effects to be 0(a) in the 
unimproved case and O(a^) in the improved one. Nevertheless we have performed 
fits on both datasets with the two Ansatze 



''The typical statistics accumulated for small lattices is of several hundred configurations. For 
the largest lattices the number of configurations ranges from around 60 at the weakest couplings to 
around 200 at the strongest ones. 





Sp(n, a/L) = (Tp{u) + p{u){a/L)'^ . 
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Table 2: Results for the step scaling function Ep, improved case. 



9 





I, 

a 






^p(5o,§) 


(50, f ) 


Sp (n, f ) 


7.7170 


6 


1 6950f26) 


133517(8) 


7527('6) 


6997(4) 


9296(9) 

vy«tyiJtyvyvty / 


7.9741 


8 


1 6950f28) 


133179f5) 


7452('6) 


6934(11) 


9305(17) 

vy«tyvyvyvy\-'-i ( 


8.3218 


12 


1 6950f79) 


132756('4) 


7353 ('41 


6858(14) 

W«WLytyt_?\ I 


9327(20) 

\j • f_y vy ^ 1 V ^ vy i 


8.5479 


16 


1.6950(90) 


0.132485(3) 


0.7266(12) 


0.6792(16) 


0.9348(27) 


7.4082 


6 


1 8811f22l 


133961 f8) 


7345 f 7) 


6773(5) 

v_/«vyi 1 fj \\j 1 


9221(11) 


7.6547 


8 


1 8811f28) 


133632(6) 


7259('7) 


6712(12) 


9246(19) 

vy«ty^^vyy -Lty / 


7.9993 


12 


1 8811f38) 


133159f4) 


0.7174(4) 


6630(13) 


9242(19) 


8.2415 


16 


1.8811(99) 


0.132847(3) 


0.7132(16) 


0.6578(14) 


0.9223(29) 


7.1214 


6 


2 1000f39) 


134423f9) 


7149('7) 


6512(5) 


9109(11) 


7.3632 


8 


2 1000f45) 


134088('6^ 


7069('6) 


6452(13) 


9127(20) 


7.6985 


12 


2 lOOOfSO) 


133599f4) 


6976f4l 


6370(14) 

vy • \y vy i vy \ -l a / 


9131(21) 

vy • ty vy \ J- 1 


7.9560 


16 


2.100(11) 


0.133229(3) 


0.6904(12) 


0.6348(12) 


0.9195(24) 


6.7807 


6 


2 4484f37) 


134994fll) 


6874('8) 


6112(5) 


8891(13) 

vy • vj vy ty -L \ vy / 


7.0197 


8 


2 4484(451 


134639(7) 


6796f7) 


6079(14) 

vy ■ vy vy i *y i -l j- i 


8945(23) 

vy • vy ty vy \ ^ vy / 


7.3551 


12 


2 4484f80) 


134141f5~) 


6711f5) 


5978(15) 

vy ■ vy ly i vy i -l vy i 


8908(23) 

vy • vy ty vy vy \ ^ vy i 


7.6101 


16 


2.448(17) 


0.133729(4) 


0.6664(12) 


0.5996(13) 


0.8998(25) 


6.5512 


6 


2.770(7) 


135327fl2) 


6628(8) 

vy • vy vy ^ vy \ w / 


5775(4) 

vy • vy 1 1 '.y \ / 


8713(12) 

vy • vy 1 J- vy \ ^ i 


6.7860 


8 


2.770(7) 


135056f8) 


6551('8) 


5753(14) 

vy • vy 1 vyvyi -l j- i 


8782(24) 

vy • vy 1 vy ^ \ ^ i 


7.1190 


12 


2 770('11'1 


134513l'5) 


6487(5) 


5704(10) 

VJ • VJ 1 Vy j: \J.\JJ 


8793(17) 


7.3686 


16 


2.770(14) 


0.134114(3) 


0.6452(14) 


0.5672(15) 


0.8791(30) 


6.3665 


6 


3.111(4) 


0.135488(6) 


0.6395(9) 


0.5427(13) 


0.8486(24) 


6.6100 


8 


3.111(6) 


0.135339(3) 


0.6356(8) 


0.5466(15) 


0.8600(26) 


6.9322 


12 


3.111(12) 


0.134855(3) 


0.6290(12) 


0.5363(15) 


0.8526(29) 


7.1911 


16 


3.111(16) 


0.134411(3) 


0.6286(9) 


0.5438(16) 


0.8651(28) 


6.2204 


6 


3.480(8) 


0.135470(15) 


0.6179(4) 


0.5058(12) 


0.8186(20) 


6.4527 


8 


3.480(14) 


0.135543(9) 


0.6129(5) 


0.5085(17) 


0.8297(29) 


6.7750 


12 


3.480(39) 


0.135121(5) 


0.6092(10) 


0.5102(15) 


0.8375(28) 


7.0203 


16 


3.480(21) 


0.134707(4) 


0.6050(10) 


0.5056(17) 


0.8357(31) 



Table 2: (continued) 
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Table 3: Results for the step scaling function Sp, unimproved case. 
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Table 3: (continued) 
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Another issue raised in ref. [2] is the number of data points which should be included 
in each fit. In that work the L/a = 6 results were dropped from the fits, being too 
far from the continuum limit. We have performed fits with all data (4-point fits) 
and also without the L/a = 6 data (3-point fits). This means that we have applied 
a total of four fitting procedures (the two Ansatze of eqs. (|5.1|5.2j) . each for a 3- and 
a 4-point fit). 

The results of these fitting procedures can be summarised as follows: 

(i) In all cases, the statistical accuracy of our result for dp is better than 1%. The 
results for the linear or quadratic coefficients p have large statistical uncer- 
tainties (up to 100%), reflecting an overall weak cutoff dependence of Ep. 



(ii) For any given lattice regularisation (i.e. improved or unimproved) and with 
any given fitting Ansatz (i.e. linear or quadratic in (a/L)), the results for cTp 
obtained by a 3-point fit are compatible to those obtained by a 4-point fit (at 
fixed coupling u). Naturally, the former have a larger error. 

(iii) For either lattice regularisation (i.e. improved or unimproved) and with any 
given number of fitting points (i.e. 3-point fit or 4-point fit) the results for 
(7p obtained by a linear fit in {a/L) are compatible to those obtained by a 
quadratic fit in {a/L) (at fixed coupling u). There is just one exception for the 
improved data at the strongest coupling u = 3.480 with a 4-point fit (agreement 
is within l.ha). The results from the quadratic fit are more accurate, due to 
the fact that the extrapolation from the range of simulated data points to the 
continuum limit is shorter in (a/L)^ than in {a/L). 

(iv) The goodness of fit is always satisfactory (x^/d.o.f. < 3) at weak and in- 
termediate couplings {u £ [0.8873,1.8811]). In a limited number of cases at 
stronger couplings the value tends to rise considerably, but this apparently 
does not depend systematically on the number of fitted points and choice of 
fitting Ansatz. In any case, given the small number of fitted data points, 
X^/d.o.f. is a goodness-of-fit criterion of relatively limited value. Instead, the 
total x^/d.o.f. varies between 1 and 2, indicating satisfactory overall quality 
of the fits. 

We conservatively consider our 3-point flt results to be our best (i.e. we drop the 
data computed at the largest lattice spacing) and opt for linear fits in (a/L) with 
the unimproved case and quadratic ones with the improved one. The results for 
these options are shown in Fig. [2j 

One could attempt to enrich this analysis along the lines of ref. we recall 
that the discretisation effects known from perturbation theory (see eq. (|4.2|) and the 
related discussion) can be divided out of the lattice SSF, by defining the quantity 

1 + n Ki[a/L) 
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The continuum limit of Sp is trivially the same as that of Sp, but the former 
quantity may approach it faster, as it has discretisation errors which are of order u^. 
However we have seen in the previous Section that ki{a/L) is always numerically 
very small. Thus the denominator of eq. H5.3|) has an imperceptible impact on Sp. 

The continuum extrapolations of Sp (obtained with improved and unimproved 
lattice actions) give results which are fully compatible both in the weak and strong 
coupling regions. At intermediate couplings we only have agreement within 1.5 a\ see 
Fig. 121 The previous fitting analysis strongly suggests that this small discrepancy, 
rather than signalling a lack of continuum limit universality, is to be attributed to 
discretisation effects not being fully under control. 

We will now corroborate this conclusion, by fitting our best results for the 
continuum SSF crp(n) with the polynomial 

N 

a^{u) = l + Y,Snu'' . (5.4) 

n=l 

In all cases the first order coefficient is fixed to its PT value, si = — 81n(2)/(47r)^. 
One-parameter fits with N = 2 yield 

S2 = -0.0029(2) (xVd-o.f. ~ 1.0) improved case , 

S2 = —0.0028(3) (x^/d.o.f. ^ 1.1) unimproved case , (5-5) 

which are not too far from the PT value S2 = —0.002031(4) of ref. [l^j. One- 
parameter fits with S2 fixed by PT and = 3 yield 

S3 = —0.00031(5) (x^/d.o.f. ~ 0.7) improved case , 

S3 = —0.00025(11) (x^/d.o.f. ^ 1.1) unimproved case . (5.6) 

The above results are compatible for the two lattice actions and thus supportive 
of universality. This analysis becomes unstable once we push it to two- or more- 
parameter fits. For instance, an = 3 fit with two fitting parameters (s2 and S3) 
yields results with errors that range between 50% and 100%, while the A = 4 fits 
with either two (53,84) or three fitting parameters (52,53,54) estimate them with 
100% uncertainty. 

Having mustered adequate numerical support for universality, we follow ref. [S] 
and calculate <Jp{u) (at fixed coupling u) by combined extrapolation of the Sp(ti, a/L) 
data from both actions, constrained to a unique continuum limit. The improved 
(unimproved) case is assumed to depend quadratically (linearly) on a/L. Results 
for dp obtained with 3- and 4-point fits are fully compatible at all couplings, while 
those for p{u) are ill-determined, as they carry up to 100% uncertainties. The good- 
ness of fit is mostly x^/d.o.f. ~ 1, except for a couple of cases where it is around 4; 
anyway its average for all couplings drops below 1. 
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The 3-point fit results for ap{u) are subsequently fitted according to eq. (|5.4j) : 
with si given by PT, the case corresponding to eq. ()5.5() gives 

52 = -0.0028(1) (xVd.o.f. ~ 1.1) combined case , (5.7) 
while that of eq. 1)5. 6|) gives 

53 = -0.00030(5) (xVd.o.f. ~ 0.8) combined case . (5.8) 

We take the results of eqs. and ((EHl) to be our best fits. In Fig. El we compare 
the LO and NLO predictions for the SSF with our discrete non-perturbative data 
and the best-fit result. 



5.2 RG running of the quark mass 

Using the functional form for fip we can compute the ratio of renormalised quark 
masses between the minimum and maximum renormalisation scales covered by our 
simulations. In order to be consistent with the notation of ref. p|, we denote the 
former by (2Lmax)~^- The ratio in question is then obtained in two steps: 
First the SSF of the gauge coupling 

^(^) = 9'{im\.,^,/L)=u ' (5-9) 

computed in ^H], is used in order to determine the correspondence between renor- 
malised couplings and renormalisation scales. This is done through the recursion 

ui = a{ui+i) , (5.10) 

with uq = ^^(l/Lmax) = 3.48 the initial value^. We note in passing that this 
procedure is based on obtaining the SSF by fitting the results of refs. ^ |2I by a 
polynomial 

N 

o-(n) ='u[l + ^a„n"j . (5.11) 

n=l 

In the present analysis we have used the = 4 series, with ai, a2 fixed from PT 
and (73, £74 resulting from the fit. 

Second the functional form for the SSF o"p is used for this sequence of couplings 
in order to compute the mass ratio from the product (cf. eq. (|3.4j) ) 

m(l/2Linax) T-f r , X1-1 

m(l/2-^^iL^a.) = n . (5.12) 



^This initial value uo = 3.48 corresponds to Lmax/^o = 0.738(16); the initial calculation was 
performed in ref. |19| while the above result is quoted in the more recent ref. |2U|. 
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Ref. 



This work 
This work 
This work 



Method 



Improved 
Improved 
Unimproved 
Combined 



M 

m(l/2Lmax) 



1.157(12) 
1.154(9) 
1.160(13) 
1.155(9) 



Table 4: Ratio of the RGI quark mass to the renormahsed quark mass at scale 
/i = l/2Lmax- "Method" refers to the procedure used in the computation of dp. 



In practice the range of scales covered by our simulations is spanned in = 7 
iteration steps. 

The final step in our calculation is the computation of the ratio of the RGI quark 
mass M to its scale dependent counterpart fn{fj,); in the quenched approximation 
this is given by |2j 



M 



22 _2, ; 

9 W 



{4tt 



-4/11 r .gif,) 

exp < - / dg 







r{g) 8 



In practice we compute the product of two ratios: 

M ( m(l/2L^ax) M 



(5.13) 



(5.14) 



m(l/2Lmax) V^(l/2-^+^W)y m(l/2-'=+iLrnax) ■ 

The first ratio on the r.h.s. is known from eq. (|5.12|) . The second ratio, which refers 
to a perturbative scale ji = l/2~^^^ L^g,^, is calculated from eq. (|5.13|) with the NLO 
perturbative values of j3{g) and T{g). 

Having described the method, we gather the relevant results (and that of ref. [2]) 
in Tabled The errors have been computed as outlined in Appendix B of ref. P]. 
The following comments are in place: 

(i) The quoted results have been obtained from the best SSF fits of eqs. (|5.6|) and 
1)5. 8() . Several other fits, such as those described in the previous subsection, 
have been tried out. In all cases the final result M/m(l/2Lmax) fluctuated 
within the quoted error, which in turn only increased slightly with increasing 
number of fitted parameters. 

(ii) Our improved result is compatible with that of ref. 0. The error is now 
smaller, due to improved statistics for the raw data on Zp. 

(iii) Compared to our unimproved result, the improved one has a smaller error. 
Recalling that statistics are roughly the same, this reflects a better system- 
atic control of discretisation effects, such as stable quadratic extrapolations in 
(a/L). 
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(iv) The result of the combined case is identical to the improved one. 



Our final result is 

M 



m(l/2Li, 



1.155(9). (5.15) 



The quoted error does not include the effect of the uncertainty in the determination 
of Ljaax/fQ, reported in ref. pO . Its contribution being roughly equal to the above 
error implies that there is no point in increasing the precision of our result unless 
the uncertainty in L^ax/fo is also reduced. 



6 Conclusions 

We have performed a very detailed computation of the step scaling function of the 
quark mass in quenched lattice QCD, employing two variants of the lattice regu- 
larisation, namely unimproved and Clover-improved Wilson actions. In both cases 
the SSF has been computed at many renormalised gauge couplings (correspond- 
ing to a wide range of renormalisation scales) and for several lattice resolutions. 
Upon extrapolation to the continuum, the SSF has turned out to be independent 
of the specifics of the lattice regularisation, providing convincing evidence for the 
universality of the continuum limit. 

The uniqueness of the continuum SSF has subsequently been used as a con- 
straint, giving us an extra handle for the control of the sensitive extrapolations to 
zero lattice spacing. The final outcome of this detailed analysis, applied to high 
statistics data, is a very precise value of the ratio M/fn{p). Far from being an aca- 
demic exercise, the increased accuracy of our result, compared to j2|, is of practical 
relevance. For example, in the context of the non-perturbative matching of Heavy 
Quark Effective Theory and QCD in finite volume, recently proposed and applied 
in refs. ||21I| . precise numerical knowledge of the functional dependence of QCD ob- 
servables on the renormalisation group invariant quark mass is of great importance. 

The analysis described in this work is currently being applied to the SSF of 
other phenomenologically interesting quantities. Preliminary results on the SSF of 
the tensor bilinear operator tpa^k'^ (relevant e.g. to some semileptonic i?-meson 
decays) have appeared in ref. [22j, while the first results on the SSF of four-fermion 
operators (related to neutral meson oscillations, Kaon decays etc.) can be found in 
ref. [2S1. 
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Figure 2: Continuum extrapolations of Sp at fixed renormalised coupling u for the 
improved action (full symbols, solid line) and the unimproved action (open symbols, 
dashed line). The L/a = 6 data points have not been included in the fits. The value 
of u increases from top to bottom and from left to right. 
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Figure 3: The step scaling function o"p(u) (full points) in the improved (top), 
unimproved (middle) and combined (bottom) cases. Shown are also the expressions 
for the step scaling function in LO (dotted line) and NLO (dashed line) perturbation 
theory, as well as our best fit to the numerical data (solid line). 
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